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Abstract 

A quantum system at equilibrium is represented by a corresponding classical system, chosen 
to reproduce the thermodynamic and structural properties. The objective is to develop a means 
for exploiting strong coupling classical methods (e.g., MD, integral equations, DFT) to describe 
quantum systems. The classical system has an effective temperature, local chemical potential, and 
pair interaction that are defined by requiring equivalence of the grand potential and its functional 
derivatives with respect to the external and pair potentials for the classical and quantum systems. 
Practical inversion of this mapping for the classical properties is effected via the hypernetted chain 
approximation, leading to representations as functionals of the quantum pair correlation function. 
As an illustration, the parameters of the classical system are determined approximately such that 
ideal gas and weak coupling RPA limits are preserved. 
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I. INTRODUCTION 



A simple atomic system in its liquid state is a prototypical strongly coupled system with 
no small expansion parameters available to simplify calculations. Under conditions where 
classical mechanics is applicable there are a number of theoretical and computational meth- 
ods available to address such strong coupling conditions. These include integral equation 
and density functional methods, and numerical techniques such as molecular dynamics (MD) 
and Monte Carlo (MC) simulation. For quantum systems there are complementary many- 
body methods to address the new features of diffraction and exchange or degeneracy effects. 
Conditions of both quantum degeneracy and strong coupling pose challenging problems with 
few choices for accurate computation of properties sensitive to both features (e.g., diffusion 
and path integral Monte Carlo (PIMC) simulation). 

One approach to describe such quantum systems is to extend the classical methods in a 



phenomenological way. For example, 
include short range diffraction effects 



VID has been used with effective potentials modified to 
l|. A related approach, wave packet MD, uses a similar 
variant of classical mechanics to define trajectories guiding wave-packets. An interesting new 
approach is to use classical integral equations for pair correlations modified to include an 
effective potential and an effective temperature Quantum effects are imbedded in these 
effective properties, while strong correlations are generated by the classical form of the 
integral equation. A surprisingly wide class of quantum systems have been described in this 
way with considerable success j3|. 

The objective here is to formalize the exploitation of classical methods by defining a 
classical statistical mechanics whose equilibrium properties are the same as those for an 
underlying quantum system. The exact definition is given in the next section, providing the 
basis for the introduction of practical approximations. One simple realization of electron 
gas as a classical system is described in the following sections. The emphasis here is on the 
definition of the mapping and its approximate evaluation. Applications will be described 
elsewhere. 
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II. DEFINITION OF THE CLASSICAL SYSTEM 



Consider for simplicity a one component system at equilibrium represented in the grand 
canonical ensemble. The grand potential | /i, 0) is proportional to the pressure and is 
given by 

I = -p{P I f,,<f>)V = -rMn J]rr^re-^(^+*-^*^W"(^)). (1) 

N 

Here Trjv denotes a trace over properly symmetrized N particle states. The thermodynamics 
for the system is specified as functions of the inverse temperature (3 and functionals of the 
local chemical potential /x(r) = fi — (p^xt (i") (an external potential has been included for 
generality) and pair potential 0. The Hamiltonian H = K + ^ contains the kinetic energy 
K and sum over pair potentials: ^ — ^ X^j^ {Qij)- 

A corresponding classical system is considered with Hamiltonian He in the same volume V 
at equilibrium described by the classical grand canonical ensemble, with inverse temperature 
f3c and local chemical potential /Xc(r) — (t>c,ext (r) . The Hamiltonian has the same form 
except the potential energy functions 0c (%) and (pcext (Qi) are different. The classical grand 
potential is defined in terms of these quantities by 

^c{Pc I l^cAc) = -pM I l^c'PejV = -lln J]^^ J rfqi..rfq^e-'^=(*-/'^'-''=('^)"W). 

(2) 

Here, Ac = {2T:^Ji? /mf'^'^ is the thermal de Broglie wavelength associated with the clas- 
sical temperature. The integration for the partition function is taken over the N particle 
configuration space. 

The classical system has undefined ingredients: the effective inverse temperature, /3c, the 
local chemical potential, /Uc(r), and the pair potential for interaction among the particles, 
0c A correspondence between the classical and quantum systems is defined by express- 
ing these quantities as functions or functionals of ^, //(r), and {qij). This is accomplished 
by requiring the numerical equivalence of two independent thermodynamic properties and 
one structural property. The first two are chosen to be the equivalence of the grand potential 
and its first functional derivative with respect to the local chemical potential. 

n (Q \ A \ - n(Q I X\ 5^c{Pc I /ic,0c) I _ 5n{B I 11,(1)) 



3 



An equivalent form for these conditions can be given in terms of the pressure and density 



Pc{/3c I /ic, 0c) = Pi(3 I /i, 0), ^c(r;/3c I l^c, 0c) = n{r;(3 | /i, 0). (4) 

These two relations provide two independent relations between Pc, /^c(r) and the physical 
variables /3 and /i(r). It remains to have a structural equivalence to relate the pair potential 
(pciQij) to 0(%), which are two particle functions. This is accomplished by equating the 
functional derivatives of the grand potentials with respect to these pair potentials 

6Qc{,/3c I /ic, 0c) I _ SQ{f3 I /i,0) , 
<50,(r,rO 50(r,rO '^''^ " 

An equivalent form for this definition is the equivalence of pair correlation functions 

9c{r, r'; /3c \ He, (pc) = 9{r, r'; /3 | /i, 0). (6) 

In this way the three equations of @ and (|6]) determine, formally, the classical parameters 
f3c, /ic, and 0c (g) as functions of (3, and functional of /i(r), and (g): /3c = /3c(/3 I 0), A^c = 
Hcij^P I /i, 0), 0c = 0c (r, I"'; /3 I At, 0)- This completes the definition of the thermodynamics 
of a classical system that is representative of that for a given quantum system. It is also 
possible to deduce from this a classical density functional theory that is representative of 
the quantum density functional theory. This will be discussed elsewhere. 



III. INVERSION OF THE MAP 



Exact representations of the classical parameters as functions/functionals of the quantum 
properties follow from classical density functional theory which has the form 

/3c/ic(r) = 2:(r I ncS-c), /3c0c (r, r') = J' (r, r' | nc, fi-c) • (7) 

Here X(r | -, ■) and J [r \ -, ■) denote formal functionals of the local density and pair corre- 
lation function. But according to the definitions ^ and ([6]) these can be replaced by their 
corresponding quantum properties to give the desired result 

Pciicij) =X{r\n,g), /3c0c (r, r) = (r, r' | n, ^f) . (8) 

This provides the formal expressions for the classical local chemical potential and classical 
pair potential in terms of properties of the quantum system, suitable for the introduction of 
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FIG. 1: (a) Dimensionless effective potential at various degeneracies (z values), (b) Effective tem- 
perature as a function of z; Asymptotic limit (dotted) implies is finite at T = 0. 

approximations. For example a good approximation to J {r \ n, g) for Coulomb systems is 
the hypernetted chain (HNC) form J' (r, r' | ri, (7) — > — ln{g (r, r')) + g (r, r') — 1 — c (r, r'), 
where c (r, r') is the direct correlation function, defined in terms of g (r, r') by the Ornstein- 
Zernicke equation. 

It remains to find an equivalent expression for Pc- This is obtained from the first definition 
of @ and the classical virial equation for the pressure, which has the form 

^ = = ^ t^'^'"^' 3c\ ^ WciJ^c, Mc, g] ■ (9) 

In the last step gc has been replaced by g according to (Q. Since /3cA*c and /3c0c are known 
from (IH]), (3c/ P is given entirely in terms of properties of the quantum system. 



IV. EXAMPLE - UNIFORM IDEAL FERMI GAS 

An interesting first illustration of this correspondence of classical and quantum system is 
the uniform ideal Fermi gas. The quantum description of this non-interacting system with 
only exchange effects is of course simple to analyze. However, the corresponding classical 
system is a true many-body problem since those same effects translate into a pair potential 
among all particles [^]. The full many-body problem must be solved to find that effective 
potential. Here that problem is addressed in terms of the approximate HNC integral equa- 
tion. Only the results will be given. In Figure (la) the classical pair potential is shown 
as a function of r/X where A is the thermal de Broglie wavelength defined in terms of the 
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FIG. 2: Classical potential times r/A showing the regularization of Coulomb potential at the origin 
and comparison with the Kelbg potential. 

quantum /3. A family of curves is given for different values of the degeneracy parameter 
z = exp for the quantum system. Of interest is the development of long range 
behavior as the temperature is lowered (large z). For this reason the mapping to a classical 
system requires a uniform neutralizing background to cancel this divergent behavior. 

The classical temperature is determined from the virial equation, with g{r) defined as 
above from the HNC integral equation. This is shown in Figure (lb). As expected for 
z < l,(3c//3 ~ 1 while approaches a finite value for T = 0. It is interesting to note that 
while the pressures of the classical and quantum systems are the same by definition, the 
internal energies are the same as well. However, this is true only if the internal energy is 
defined thermodynamically and not mechanically. 

V. ELECTRON GAS 

A practical classical representation for the uniform interacting electron gas can be ob- 
tained in a similar way. The classical potential is divided into its ideal gas contribution 
iPc'PcY^^ plus a contribution from the Coulomb interactions A: I3c4>c = (/^c^c)'"'^^ + A . An 
approximation for the latter is obtained by requiring the correct RPA form at weak cou- 
pling. In the classical context this means that the direct correlation function is given by the 
classical effective potential. This is obtained from the Ornstein-Zernicke equation with the 
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classical gdr) replaced by the quantum g{r), according to the definition (|6]), evaluated in 
the RPA. The final result for the effective classical potential is 



(10) 



Here S^^\k) is the ideal gas static structure factor and S^^'^{k) is the RPA structure factor. 
It is possible to show that the familiar Kelbg potential [l| is recovered in the weak coupling, 
low density limit. Figure (2) illustrates the quantum corrections to the bare Coulomb form 
for = 5, and t = T/Tp = 5 (where Vg is the usual ion sphere radius relative to the Bohr 
radius, and Tp is the Fermi temperature). 
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